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04 ■ Abstract: We discuss recent progress in various problems related to bound-state quantum 

f—j ' electrodynamics: the bound-electron g factor, two-loop self-energy corrections and the laser- 

' dressed Lamb shift. The progress relies on various advances in the bound-state formalism, 

including ideas inspired by effective field theories such as Nonrelativistic Quantum 
Electrodynamics. Radiative corrections in dynamical processes represent a promising field for 
further investigations. 
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1. Introduction 

o ■ 

l/^ . This is a brief summary of a number of recent advances in our understanding of bound-state quan- 

' tum electrodynamic (QED) effects. The topics are (i) two-loop corrections to the bound-electron g 

fH , factor, (ii) higher-order two-loop corrections to the self-energy of a bound electron, and (iii) the laser- 

Q^' dressed Lamb shift. 

^ • The first two of these rather diverse topics are related to two-loop effects. The investigation of these 

(U [ is simplified considerably by the use of effective field-theory techniques inspired by Nonrelativistic 

• QED (NRQED) [1-3]. The Wilson coefficients multiplying the effective operators in the NRQED 

Lagrangian are matched against those of the full relativistic theory, providing a simplified framework 
for the calculation of bound-state effects. Scale-separation parameters such as the photon mass /i are 
cancelled at the end of the calculation. The analysis of higher-order corrections to the g factor of 
' the bound electron is simplified further by a transformation to the length gauge, which results in a 

lesser number of terms to be considered than would be necessary in the velocity gauge. This fact has 
inspired the development of Long-wavelength QED (LWQED) [4], a theory which is obtained after 
Power-Zienau and Foldy-Wouthuysen transformations of the first-quantized Lagrangian; the second 
quantization is carried out by formulating the path integral. Consequently, an improved understanding 
and a tremendous simplification results for the calculation of a number of QED corrections for bound 
states, such as the g factor and higher-order corrections to the self-energy. 

A further field of recent studies has been concerned with the interaction of a laser-dressed bound 
electron with the radiation field [5-8]. This process entails corrections which can only be understood 
if the analysis is carried out right from the start in the framework of the laser-dressed states, which are 
the eigenstates of the quantized atom-laser Hamiltonian in the rotating-wave approximation [9]. 
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2. Bound-electron g factor 

In this section, we briefly summarize the results of a recent investigation 1 10| of the bound-electron 
g factor, which is based on NRQED. The central result of this investigation is the following semiana- 
lytic expansion in powers of Za and In(Za) for the bound-electron g factor (nS state) in the non-recoil 
linait, which is the limit of an infinite nuclear mass: 

^ V ' 

Breit (1928), Dirac theory 

^ V ' 

one-loop correction 

+ (^)' {-0.656958 (l + ^-^^ + ln[(Zo)-] + + OiZaf] 

^ V ' 

two-loop correction 



+ 0{a^). (1) 

This expansion is valid through the order of two loops (terms of order are neglected). The 
notation is in part inspired by the usual conventions for Lamb-shift coefficients: the (lower case) 
a terms denote the one-loop effects, with auj denoting the coefficient of a term proportional to 
a (Za)^ In-* [{Za)~'^]. The b terms denote the two-loop corrections, with hkj multiplying a term pro- 
portional to {Za)'' In-' [{Za)~'^]. In [10], complete results are derived for the coefficients a4i, 040 
and 641. 

In general, the expression corresponding to (1) for a free electron is obtained by letting the param- 
eter Za ^ in every term of the loop expansion (expansion in powers of a). In this limit, the known 
free-electron two-loop result is recovered [11-17]. 

Up to the relative order {Za)"^, the free-electron contribution in one-, two-, and higher-loop order 
is multipUed by a relative factor 

1 + (Zaf 020 = 1 + {Zaf 620 = 1 + ■ (2) 

This result consequently holds for the three-loop and the four-loop term not shown in Eq. (1). The 
appUcability of the relative factor (2) to the two-loop term, valid through (Za)"^, had been stressed 
previously in [18]. The result in Eq. (2) had been obtained originally in [19-23] (for the IS state). As 
is evident from Eq. (1), the correction of relative order {Za)^ is different on the level of the tree-level 
diagrams and reads g{nS) ~ 2 x [1 - (Za)^ /{3n'^)]. 

Explicit results for the coefficients in (1), restricted to the one-loop self-energy, read [10] 

32 

a4i{nS) = —, (3a) 
^40 (nS) = -- — -- In A;o (nS) --Inks {nS) . (3b) 
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Fig. 1. Feynman diagrams for the two-loop self-energy corrections to the 
bound-electron g-factor. 



Here, lnfco(n5) is the Bethe logarithm for an nS state, and lnks{nS) is a generalization of the 
Bethe logarithm to a perturbative potential of the form Vacuum polarization adds a further n- 
independent contribution of (—16/15) to 040 [24]. The Bethe logarithms for IS and 2S [25] read 



lnA;o(lS') = 2.984128 555, 
lnfco(2S') = 2.811 769 893, 

and the corresponding values for In ks read [10] 

lnfc3(15) = 3.272 806 545, 
lnA;3(25) = 3.546 018 666. 

The quantity In ^3(718) is defined as, 



Jo 



dk k^ 



E-Hs-k E-Hs-k 



r\4>) 



In 



H 5 
(Z^^ 6 



Inks 



(4a) 
(4b) 

(5a) 
(5b) 



(6) 



where the ultraviolet cutoff e is to be understood in the sense of [2], and the matching of the nonco- 
variant cutoff e to the covariant photon mass fj, is given as (see [26], pp. 361-362) 



In 



-) 

m J 



In 



(-) 
Vm/ 



(7) 



However, this replacement is not unique and the constant term depends on the actual form of the 
integrand. A different replacement has to be used for some of the low-energy photon corrections to the 
g factor [10]. 
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The results for the two-loop coefficients read 
biiinS) = —, (8a) 

64o(lS) = -18.5(5.5) . (8b) 

Here, the result for 640 is an estimate based on an explicit calculation of a large contribution due 
to low-energy virtual photons, and an estimate of the remaining, unknown contribution due to high- 
energy virtual photons. The dominant logarithmic two-loop term 641 is caused exclusively by the two- 
loop self-energy diagrams in Fig. 1 alone. The other two-loop diagrams, which include closed fermion 
loops, can be found in Fig. 21 of [27]. The logarithmic term 641 is, however, exclusively related to the 
gauge-invariant subset displayed in Fig. 1. 

The newly calculated 041, 040 and 641 are bound-state corrections to the electron gi-factor of order 
(a/Tr) (Za)'^ and (a/ir)'^ (Za)^, multiplied by logarithmic terms. These corrections are (at Z = 1) for- 
mally of order and and therefore of the same order of magnitude as the tenth- and twelfth-order 
corrections to the free-electron anomaly, which barely are of experimental or theoretical significance 
at the current level of accuracy. One may therefore ask why these binding corrections are of any phe- 
nomenological significance. The reason is that at somewhat higher Z, the situation changes drastically, 
due to scaling of the binding corrections. In addition, due to numerically large coefficients and 
logarithmic factors, the "hierarchy" of the corrections changes drastically. Roughly, one may say that 
at Z = 1, the bound-electron anomalous magnetic moment is approximately independent of binding 
corrections of order (a/w) {Zaf' and higher, whereas for higher Z, the situation is reversed, and the 
binding corrections to the one- and two-loop contributions are numerically much more significant than 
the higher-loop free-electron corrections. This "transition from free to bound-state quantum electrody- 
namics" as a function of Z is a somewhat pecuUar feature of the bound-electron g^-factor. 

For, example, we consider the ratio 

n (Z) = ^ ——^ , (9) 

(a/7r) 

which gives an order-of-magnitude estimate for the the ratio of the one-loop self-energy binding cor- 
rection to the eighth-order anomalous magnetic moment of the free electron. We have 

ri(Z=l)«2, ri(Z = 10)«104>l. (10) 

For the two-loop logarithmic binding correction, we have 

{Zaf ln[(Za)- 
(a/Tr)^ 
and consequently 

r2(2 = 1) « 0.2 < 1, r2(2 = 10) « 2 X 10^ > 1. (12) 

As is evident from these considerations, the [Za)'^ one-loop and two-loop binding corrections are 
roughly of the same order of magnitude as the highly problematic four-loop corrections [28-30] for 
the free electron. However, the situation changes drastically even at very moderate nuclear charge 
numbers, and the binding corrections to the one-loop and two-loop contributions become dominant 
over the higher-loop effects. 



T^{Z) = ^ (11) 
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The NRQED one-loop calculation [10] is divided into three parts, the first of which entails fully 
relativistic form-factor corrections including lower-order terms, the second of which corresponds to 
a spin- dependent scattering amplitude, and the third of which is a low-energy Bethe-logarithm type 
correction that contains In fcg and In k^. The one-loop correction in Eq. (1) can therefore be written in 
a natural way as Sg^^^ = g[^^ + g^^ + g^^\ where 



(1) a 

TT 



, (Za)2 (Za)4 /7 5 16 , 



6 24n 



(1) a (ZaY f , 16 , \ 
IT n-^ \ 9 / 



(1) _ a {ZaY 32 



^3 = - 



fi 5 In/co 3 



(13a) 

(13b) 
(13c) 



The new contribution of order (Za)'^ can be compared with the numerical results for the self- 
energy correction [31, 32] complete to all orders in Za. Assuming correctness of the logarithmic term 



in Eq. (13), a fit to numerical data yields (IS) = -10.2(1) and a^o^(2S) = -10.6(1.2) forthecon- 
stant term, in excellent agreement with the analytic results which read (IS) = — 10.236 524 318(1) 
and 4o^(2S) = -10.707 715 607(1) . 

Having verified the consistency of the analytic [10] and numerical results [31, 32], an interpolation 
procedure [33] may now be used to extract a more accurate theoretical prediction at low and interme- 
diate nuclear charge numbers, if combined with numerical results at higher Z [31]. Thus, the results in 
Eqs. (3) and (8) may be used in order to infer improved theoretical predictions for the bound-electron 
g factor, notably in the experimentally important special cases of hydrogenlike carbon |34| and oxy- 
gen [35]. Alternatively, the improved status of the theory may be used in order to infer a more accurate 
value of the electron mass. Specifically, the value from the carbon measurement [34], using the new 
theory, reads 



m(^2c5+) = 0.000 548 579 909 41 (29)(3) u . 



(14) 



The first error comes from the experiment [34], and the second error corresponds to the theoretical 
uncertainty. The conclusion is that a further improvement of the experiment could lead to a much 
better determination of the electron mass; the new theory provides room for at least an improved 
determination by one order of magnitude. 

For the calculation of yet higher-order binding corrections to the one-loop and two-loop contribu- 
tions, a detailed understanding of the two-loop form-factors, including their slopes, is required. The 
most recent calculations of these effects, in both dimensional and photon-mass regularizations, can be 
found in [36-38]. 



3. Two-loop Bethe logarithms 

As is well known [39^6], the two-loop Lamb shift AiJ^^^ , in the limit of an infinite nuclear mass, 
may be written as 

AE^^^ = (^Y (^")>^' g(Za). (15) 
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Fig. 2. The two-loop corrections to the Lamb shift in hydrogenlike systems fall naturally into three gauge- 
invariant subsets, which are the pure two-loop self-energy terms (2LSE), the vacuum-polarization correction to 
the virtual-photon line in the self-energy (SVPE), and the self-energy vacuum-polarization and pure two-loop 
vacuum-polarization corrections (SEVP). The two-loop Bethe logarithm is a numerically large correction to the 
nonlogarithmic term of order o? {Zcxf and it is exclusively related to the 2LSE subset; however, for a complete 
result in this order, contributions from the other gauge-invariant subsets will have to be considered as well. 



For S states, the dimensionless function H{Zq.), has a semianalytic expansion of the form 

H{Za) =Bio + {Za) B50 

+ {Zaf [^63 ln3[(Za)-2] + \r^[{Za)-^] + Bgi \n[{Za)-^] + Beo] , (16) 

where we ignore higher-order terms, and upper case is used for the -B^ coefficients that multiply 
terms of order {Zaf In^ [{Za)^'^] nie c^. The coefficients, restricted to the two-photon self-energy 
diagrams (Fig. 2), read as follows 

B^f^^\nS) = -^ = -0.296296, (17a) 

B^f ^^\lS) = 1^-^ ln(2) = -0.639 669, (17b) 

B^r-\nS) = B^r-\1S) + ^ g + ^ - i - ln(n) + *(n) + . (17c) 
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The n-dependence of Bgi has been clarified in [46,47] 

Bg^^^^ (nS) = B^f (IS) + ^ [N{nS) - iV(lS)] 

where C = 0.577216 ... is Euler's constant, 'i'{n) is the logarithmic derivative of the Gamma function, 
and N{nS) is related to a correction to the Bethe logarithm induced by a Dirac-delta potential. Explicit 
values for N{nS) can be found in [47] (n — 1, . . . , 8). 
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Fig. 3. The integration regions for the two virtual photons 
in the two-loop self-energy problem comprise a low-energy 
regime with two low-energy photons, which gives rise to bL- 
The middle-energy regions with one low-energy and a one high- 
energy photon give rise to 6a/. The high-energy contribution 
&F + 6ir is as yet unknown. 



The Bqo coefficients are the sum of several contributions 

4f "^""^ (nS) = bL + hu + {bp + + &vp} , (19a) 

= 6l('t-S) ~ two-loop Bethe logarithm, two soft photons , (19b) 

bu — bM{nS) — —N{nS) ^ one soft, one hard photon , (19c) 

bp ^ soft electron momenta, two hard photons , (19d) 

bn ^ hard electron momenta, two hard photons , (19e) 

byp ~ vacuum-polarization corrections . (19f) 
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Only the terms + bM are currently known [48, 49] (see also Fig. 2). The contributions in curly 
brackets in Eq. (19a) remain to be evaluated. However, an estimate for the total value of Bqq may be 
obtained, 

Beo{nS) = 6l + &m + {&f + + J'vp} , 

Beo{nS) « 6l + &m ± 15%. (20) 

This estimate [48, 49] is based on corresponding one-loop calculations, where the low-energy virtual 
photons give the by far dominant contribution to the constant term [2]. The results for the two-loop 
Bethe logarithms of S states read [48,49] 

bLilS) = -81.4(3), (21a) 



bL{2S) = 


-66.6(3) , 


(21b) 


6l(35) = 


-63.5(6) , 


(21c) 


6l(45) = 


-61.8(8), 


(21d) 


bL{5S) = 


-60.6(8) , 


(21e) 


bues) = 


-59.8(8) . 


(21f) 



A few clarifying remark might be in order. The Bqo coefficient multiplies a correction of order 
a^{Za)^, which is effectively an order-a^ contribution to the energy levels of hydrogen (Z = 1). 
In order to complete the calculation at this order of magnitude, it would also be necessary to consider 
the four-loop Dirac form-factor slope of the electron, as well as the three-loop binding correction of 
order (Za)^ rUe (? [50]. The three-loop slope has recently been evaluated in [51], completing the 
theory of energy levels in hydrogen up to the order of aJ . 

4. Laser-dressed Lamb shift 

In the recent past, seminal advances have been obtained both in the techniques of high-precision 
spectroscopy (e.g., [52]), and in the coherent preparation and manipulation of media by external elec- 
tromagnetic fields [53, 54]. Thus it is desirable to study the bound electrons interacting simultaneously 
both with the quantized radiation field and with an external driving field. An accurate theory of such 
systems, including all dynamic effects, might eventually open a possibility for a whole new class of 
high-precision experiments, provided that technical problems related to the required highly accurate 
intensity stabilization of the laser (and others) can be solved. Traditionally, radiative and relativistic 
corrections are treated with methods of QED, whereas studies related to the dynamical nature of the 
interaction of matter with driving laser fields are the domain of Quantum Optics (QO). Obviously, 
a treatment of bound electrons in the presence of both the radiation field and external driving fields 
requires a combination of ideas from both subject areas: While, a priori, the essential-state approx- 
imation of QO [53] is not sufficient to obtain the accuracy of QED, a perturbative treatment of the 
interaction of the bound electron with a strong external (laser) field as in QED is hopeless because of 
the large coupling parameter. 

In [5-7], an atom with two relevant energy levels driven by a strong near-resonant monochromatic 
laser field is studied as the easiest model system for the above problem. The incoherent part of the 
resonance fluorescence spectrum emitted by this system in QO is known as the MoUow spectrum, 
where the coupling strength is characterized by the the Rabi frequency f2 defined as (^ = eo = c = 1) 

0. = -q{e\x- e\\g)E\,, (22) 
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for a driving laser field Ez,{t) = Sh^L cos(wLi) with frequency wl, macroscopic classical amplitude 
£h and polarization cl- Here, q = — |g| is the elementary charge. The corresponding coupling constant 
5l for the interaction of a quantized driving laser field with the main atomic transition is defined by 

gj^ = -q{g\ej^-x\e)£i^\ (23) 

where 5^^^ = y^aJh/^V is the electric laser field per photon and V is the quantization volume. The 
matching of the electric field per photon with the corresponding classical macroscopic electric field f l 
is then given by 

2v/n+T4^^ < — > £l. (24) 

If O is larger than the natural decay width F of the transition, then the Mollow spectrum approxi- 
mately consists of one central and two sideband peaks of Lorentzian shape, which are located symmet- 
rically around the driving laser field frequency. The sideband peaks are shifted from the driving field 
frequency by the generalized Rabi frequency JIr = + A^, where A = wl — wr is the detuning of 
the laser field frequency uih with regard to the atomic transition frequency cor. The shape of the Mollow 
spectrum may easily be explained in terms of the dressed states, which are defined as the eigenstates 
of the quantum optical interaction picture Hamiltonian describing the matter-light interaction. In trans- 
ferring to the dressed state picture, the interaction with the driving laser field is accounted for to all 
orders. 

Thus, when evaluating radiative and relativistic corrections to the Mollow spectrum, it is natural to 
start the analysis from the dressed-state basis as opposed to the unperturbed atomic bare-state basis. 
In [5, 6], it was shown that this distinction in fact has to be made. It is not sufficient to modify the 
energies (which enter in the formula for the dressed states) according to the usual bare-state Lamb shift 
in order to obtain the correct result for the corrections to the MoUow spectrum. Instead, at nonvanishing 
detuning and nonvanishing Rabi frequency, a treatment starting from the dressed-state basis leads to an 
additional nontrivial correction term. This term gives rise to a shift of the MoUow sidebands relative to 
the central peak given by 



= =pC , " , (25) 



where 

C = ^ ln[(Za)-2] (26) 

is a dimensionless constant. Here, the notation (.)g and (.)e denotes the expectation value evaluated 
with the ground or excited atomic state, respectively. 

Inspired by the interpretation of the bare Lamb shift correction in terms of a "summed" shift as 
in [5, 6], this additional correction can be interpreted as a modification to the Rabi frequency: 

SZJ^f^ = ± (1 _ c)2 + A2 - + A2) , (27) 

with 5iX7j?^ « Slo^^ because of the smallness of the correction. It should be noted that this interpreta- 
tion in terms of a summation is not trivial and was shown to be valid up to first order in the correction. 

In [5,7], the leading relativistic and radiative corrections up to relative orders (Za)^ and a{Za)'^, 
respectively, have been evaluated, as well as all other relevant correction terms up to the specified 
order of approximation. It turns out that all corrections may be interpreted as either corrections to the 



20CM NRC Canada 



10 



Canadian Journal of Physics Voi. 1 , 2004 



Rabi frequency O or as corrections to the detuning A, such that one can define the fully corrected 
generaUzed Rabi frequency Cl^^ by 

= ^n^-{i+ci^y+{A-Aii!,y . (28) 

^ (i) 

Here, OJ;^^ contains all corrections to the Rabi frequency, namely the relativistic and radiative correc- 
tions to the transition dipole moment, field-configuration dependent corrections, higher-order correc- 
tions to the self-energy, and corrections to the secular approximation. A'fJ^ consists of all corrections 
to the detuning, i.e. the bare Lamb shift, Bloch-Siegert shifts, and off-resonant radiative corrections. 
The superscript (j) indicates the dependence of the result on the total angular momentum quantum 
number. 

Equation (28) summarizes the main result of this study: In the presence of driving laser fields, the 
usual bare state Lamb shift of the atomic states is augmented by additional correction terms. These in 
part depend on the laser field parameters fl and A, which span a two-dimensional parameter manifold 
determining the actual value of the dynamical Lamb shift. 

A promising candidate for the experiment are the hydrogen IS 1/2 <-» 2P1/2 and 1 6*1/2 ^ 2P3/2 
transitions. We consider here as a specific example the IS1/2 ^ '^Py2 transition with fl = 1000 • ri/2 
and A = 50 • Ti/^ as the laser field parameters. The Rabi frequency is shifted with respect to the 
leading-order expression fi^ = \/0^"+A^ by relativistic and radiative corrections as follows, 

± (n'^c^^^ - Or) = ±738.282(60) • 10^ Hz . (29) 

This driving laser field parameter set is expected to be within reach of improvements of the currently 
available Lyman-a laser sources [55] in the next few years. The corresponding result for the ISy^ 
2P3/2 transition with Q = 1000 • Ta/j, A = 50 • Ts/j is 

± (^f2[f/'' - Or) = ±734.871(60) • 10^ Hz . (30) 

All given uncertainties are due to unknown higher-order terms [7]. 

By a comparison to experimental data, one may verify the presence of dynamical leading- 
logarithmic correction to the dressed-state radiative shift in Eq. (25), which cannot be explained in 
terms of the bare Lamb shift alone. This allows to address questions related to the physical reality of 
the dressed states. On the other hand, the comparison with experimental results could also be used 
to interpret the nature of the evaluated radiative corrections in the sense of the summation formulas 
which lead to the interpretation of the shifts as arising from relativistic and radiative corrections to the 
detuning and the Rabi frequency. 

5. Conclusions 

One of the obvious conclusions to be drawn from the recent advances in bound-state quantum 
electrodynamics is as follows. A widespread opinion has been invahdated which suggested that two- 
loop corrections to the bound-electron g factor, and two-loop corrections to the Lamb shift in higher 
order might be a severe, if not insurmountable, obstacle against further theoretical progress. Quite to 
the contrary, the recent advances have shown that an understanding of these effects is feasible to a 
good accuracy. While some important contributions remain to be evaluated, the further program (e.g., 
the evaluation of the remaining high-energy parts) is clearly defined, and decisive first steps toward a 
much improved understanding have been accomplished. 

Further advances are possible both on the experimental as well as on the theoretical side: con- 
cerning the g factor, the recent theoretical progress allows for an order-of-magnitude improvement of 
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the value for the electron mass based on a potential new measurement alone. Regarding the Lamb 
shift in hydrogen and low-Z hydrogenlike systems, one has recently gained an improved understand- 
ing of the binding corrections of order iZaf" In-' \{Zol)^'^\ (? (where j may assume the values 
j — 0, . . . , 3). Recent numerical investigations in the regime of intermediate nuclear charge num- 
bers [56] have also contributed toward an improvement of our understanding. The extrapolation of the 
low-Z results by deferred Pade approximants [57] suggests a rather good general consistency of both 
approaches, while some issues regarding the consistency of the analytic and numerical approaches 
remain to be addressed [56]. 

Other progress concerns the modifications of radiative corrections in dynamical processed as op- 
posed to ^-matrix energy shifts. The laser-dressed Lamb shift [5-8 1 is a dynamical correction to the 
dressed-state [9] quasi-energies. The self-energy, in this case, gives the by far dominant effect. While 
the bulk of the laser-dressed Lamb shift can be understood in terms of a radiative correction to the de- 
tuning, which is taken into account in a natural way by evaluating the detuning in terms of the observed 
(low-intensity) resonance frequency of the transition, a few tiny shifts persist which can only be under- 
stood if the system is treated in the dressed-state picture right from the start. The dynamic Lamb shift is 
an effect which depends on two parameters that determine the dynamics of the system: (i) the detuning 
and (ii) the Rabi frequency. Therefore, the dynamical Lamb shift could be mapped out as a function of 
these parameters in a possible experiment. A rather promising candidate for a possible measurement 
would be based on the hydrogen lS-2Pi/2 transition [5]. However, a very attractive alternative would 
be provided by a forbidden Ml transition [58] in Ar XIV, 2s^ 2p ^Pi/2-^P3/2' provided the many-body 
QED effects can be treated to sufficient accuracy. The system in question is described very well by a 
two-level formalism, and the resonance line width is small. 
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